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ABSTRACT 

We show how the viscous evolution of Keplerian accretion discs can be under- 
stood in terms of simple kinetic theory. Although standard physics texts give a simple 
derivation of momentum transfer in a linear shear flow using kinetic theory, many 
authors, as detailed by Hayashi & Matsuda 2001, have had difficulties applying the 
same considerations to a circular shear flow. We show here how this may be done, and 
note that the essential ingredients are to take proper account of, flrst, isotropy locally 
, in the frame of the fluid and, second, the geometry of the mean flow. 
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INTRODUCTION 



o 
o 

^^i; Accretion discs play a central role in a wide range of astronomical environments, mediating the gas flows in the vicinity of 
O ' object as diverse as AGN, binary stars and protostars (Pringle, 1981). In an accretion disc, the predominant flow is a circular 
shear flow, with angular velocity n{R) a function of radius R from the central object. Accretion takes place because of the 
action of some form of dissipation which releases the free energy of the shear flow as heat, and so allows the disc material to 
fall deeper into the potential well of the central object. Simple physical energy arguments {e.g. Lynden-Bell & Pringle, 1974) 
indicate that the dissipative process must take the form of a stress which transports angular momentum outwards. Because 
■ the free energy of a circular shear flow is zero only if dQ,/dR — 0, it follows that the relevant element of the stress tensor must 
$_( ' be of the form 

■ 

oc ~dn/dR, (1) 

where is the azimuthal coordinate. 

This can be deduced from the standard derivation of Navier-Stokes stress to be found in the fluid dynamics textbooks. 
However, in a recent paper, Hayashi & Matsuda (2001) have drawn attention to the fact that attempts to provide a physical 
explanation of the above result in terms of simple kinetic theory have resulted in failure. 

The simple kinetic explanation given in basic physics text books for the effect of viscosity on a simple linear shear flow, of 
the form u — (0, U{x)) in Cartesian coordinates, relies on the fact that the kinetic particles conserve linear momentum between 
collisions. Thus particles crossing some flducial plane x = Xo tend to mix up and smooth out the momentum distribution of 
the fluid and so give rise to a stress of the form Tyx (x ~dU/dx. 

However, as Hayashi & Matsuda (2001) point out, it is in the application of these simple concepts to a circular shear 
flow that the problems seem to arise. The simple generalisation that, in a circular shear flow, the kinetic particles conserve 
angular momentum j — R^Q between collisions would imply, taken at face value, that the the movement of particles across 
some fiducial circle R = Ro would tend to try to mix up and smooth out the distribution of angular momentum of the fluid, 
and thus that the stress would be proportional to —dj/dR (see e.g. Madej & Paczynski 1977). From the arguments given 
above, this is clearly wrong. Not only would this predict a stress in the case when the shear dQ/dR is zero, but for a standard 
Keplerian accretion disc for which j oc R^^^ it would transport angular momentum inwards rather than outwards. As detailed 
by Hayashi & Matsuda (2001) attempts to get round this and to produce the 'correct' result have only succeeded by making 
mathematical errors in the derivation. 

From all these problems, Hayashi & Matsuda (2001) conclude that although what they call the derivation 'with mathe- 
matical rigour' (i.e. the usual Navier-Stokes argument) gives the correct answer, in order to obtain the correct answer using 
kinetic theory one must take account of such complications as Coriolis force. In this paper, we shall show that, although it 
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is obviously necessary to include Coriolis force if one works in a frame co-rotating with the flow, one can obtain the correct 
result from straightforward kinetic theory in the inertial frame. 

Before wo do so, it is instructive to return to the 'mathematical' relationship between stress and strain derived in the 
standard fluid textbooks (see, for example, Batchelor, 1967, Section 3.3). The basic point we wish to make is that the 
standard Navier-Stokes expression for momentum transfer {i.e. stress) is based on a simple physical argument. The argument 
may involve the use of tensors, which physicists tend to meet in courses on mathematical methods, but the argument itself 
is not a mathematical one. The point is that stress (momentum transfer) in a fluid (or in a solid) can be expressed as a 
second order tensor. This is a physical (and not a mathematical) statement, in exactly the same sense that the statement 
that a velocity is a first order tensor (i.e. a vector) is a physical, and not a mathematical, statement.^ For a fluid, the physical 
Ansatz is simply that the stress tensor must be physically related to the rate of strain tensor (which is a second order tensor 
which derives from the first derivatives of the velocity field, and thus incorporates information about the shear). The simplest 
assumption is that relationship between these two tensors is through a (fourth order) tensor which is isotropic. It is this 
assumption of the isotropy of the relationship, which is based on the physical assumption that the fluid itself is isotropic, 
which gives rise to the standard Navier-Stokes expression for the viscosity.^ 

It is important to realise that these tensors (i.e. scalars, vectors, second order tensors) exist as physical quantities, 
independent of any coordinate system. The statement that there is a relationship between two of them is a physical statement. 
It is only when one has to calculate a particular element of the stress tensor, for example corresponding to linear momentum 
transfer in a linear shear flow, or to the angular momentum transfer in a circular shear flow, that one has to evaluate coordinate 
dependent expressions, which can get mathematically complicated. But when one does this, one finds that the flux of linear 
momentum in a linear shear flow just depends on —dU/dx, and that the angular momentum fiux in a circular shear flow 
just depends on —dQ/dR. However, this result also enables us to draw another conclusion. The reason for the difference 
between the terms dU/dx and dfl/dR is due solely to the difference between the coordinate systems. That is, it comes from 
geometry alone. This implies that when looking for differences in derivations for simple kinetic theory formulae between the 
linear and circular shear flows, we need only concern ourselves with geometry, and not with dynamical complications such as 
Coriolis force. In addition, we also need to take note of the fact that the Navier-Stokes expression does depend critically on 
assumptions about isotropy of the fluid. Thus we should expect to have to make a similar assumption about the properties of 
our kinetic particles. 

2 KINETIC THEORY 

In this Section we compute the relationship between stress and (rate of) strain for two simple shear flows using kinetic theory. 
To keep the concepts and the algebra simple, we make a number of simplifying assumptions. We work in two dimensions only. 
That is, we assume that the kinetic particles (assumed to be identical, with mass m) move only in a two-dimensional plane. 
We assume that the net effect of the scattering processes within the fluid is that these particles are emitted at a constant 
rate at each point of the fluid, and that each particle is emitted with an identical velocity, c, relative to the local fluid. ^. We 
represent by iV(A) the number of particles emitted per unit area per unit time that travel a distance greater than A before 
colliding with another particle and we assume that N{X) is independent of position. The requirement of isotropy, as discussed 
above, can now be imposed by making the assumption that this emission takes place at each place in the fluid isotropically 
in the frame of the fluid at that place. 

2.1 Plane shear flow 

In Cartesian {x, y) coordinates we let the background fluid flow be of the form 



We consider the stress acting on a line element, length dl, centred at point S, at position (a::o,0), and lying in the 
direction, that is, with unit normal in the a:— direction. For a linear shear flow, we consider the case 



^ The mathematical complications come in when one lias to specify a particular tensor (or vector) by specifying the coordinate repre- 
sentation of that tensor in a particular coordinate frame, and get worse when one has to specify the representation of that same tensor 
in some other coordinate frame. 

This simplest assumption is not necessarily the correct one. For example it is possible that the viscosity might depend on the absolute 
orientation of the shear in some inertial frame {e.g. the stress due to convection in a rotating medium - see Kumar, Narayan & Loeb 1995, 
and references therein; the stress due to non-isotropic mixing - Bretherton & Turner, 1968; or the stress induced by the magneto-rotational 
instability - Torkelsson et. al. 2000). 

^ It is of course straightforward to generalise the analysis to take account of the particles being emitted with a distribution of velocities. 



u={0,U{x)). 



(2) 



U{x) = Uo + {x- xo)U', 



(3) 
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Figure 1. (a) Schematic of linear shear flow showing emission point, E, line element S and path of particles, ES, in frame of S. The 
dashed arrow denotes the velocity of of the mean flow at E with respect to that at S, the bold arrow the random velocity component 
and the open arrowhead denotes the particle resultant velocity in the frame of S. Note that in reality (although not in the Figure) the 
length of the bold arrow is much greater than that of the dashed arrow, (b) Schematic of linear shear flow showing the patch around E 
that is the area element for integration of the viscous stress. 



where U' = dU/dx is a constant. Note that although we consider here a linear shear flow with constant shear, any linear shear 
flow can be treated as a constant shear for a: « xo and with U' = dU/dx, evaluated at x = xo- 

For convenience, we work in the frame comoving with the mean fluid flow at S. This frame is an inertial frame, and in 
this frame the mean flow is given by 

u = (0, {x - xo)U'). (4) 

We now focus on particles that are emitted from a point E that is located at a distance A from S and where the line 
SE makes an angle a with the negative x— axis (see Figure la). The critical point is that in the frame of the fluid at E the 
particles that travel along ES are emitted at an angle 9 to the a— axis where the relationship between 6 and a is simply 
deduced by considering the ratio of the x— and y— components of the particle velocity in the rest frame of S, i.e. 

csinO — U' \cosa 

tana = . (5) 

ccos9 

We consider the limit that the shear velocity across a typical mean free path is much less than the random particle 
velocity {i.e. U'\ <^ c), so that (to 0{U'\/c)) we may write 

tan ~ tan a + U'X/c, (6) 
cos^ ~ cosa — sinacos^ a J7'A/c, (7) 
and 

sin^ ~ sina + cos^ a ?7'A/c. (8) 

In order to compute the rate of arrival of y-momentum at S due to particles originating near E, we consider a patch 
around E that is a portion of an annulus centred on S (radius A, thickness dA) where the patch subtends an angle da at S 
(see Figure lb) . The rate of emission of particles from this patch that travel far enough before colliding to be able to reach S 
is simply N{X)\d\da. Such particles emanate isotropically from E in the rest frame of the fluid at E. Thus the fraction, /, of 
such particles that impinge on the line element dl at S is 

_ dl cos a , , 

27rL ' ^ ' 

where L is the distance traveled by the particles between E and S in the rest frame of the fluid at E. The relationship between 
L and A is simply ascertained by noting that the relative velocity between E and S is zero in the a;— direction and hence that 
the distances travelled in the x— direction between E and S in the two frames is thus equal (i.e. LcosO = A cos a) from which 
we deduce, to first order in A, 

L ~ A(l — sinacosaf7'A/c)~^, (10) 
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and thus that 

d/ cos q(1 — sin Q cos q[/'A/c) , , 

Finally, the y— velocity of each particle emitted from around E is 
Uj, = — )7'A cos a + csin 6', (12) 
which (using equation |5J may be approximated as 



Uy 



' C 



sin Q — t/'A cos a sin^ a. (13) 



Thus the total rate of arrival of y-momentum at S from the patch at E is 

Py = mvyfN{\)XdXda. (14) 

Substituting for / and Vy from equations II II and 1131 we may obtain, to the order of our approximation, the total rate of 
arrival of y— momentum at S by integrating over a and A, in the form 

Py ^ / cos a(l — sin a cos Q(7'A/c)(csin Q — ?7'A cos Of sin^ Q)dQ: dA. (15) 

Jo 27r J_^^2 

Since we are only considering the momentum flux due to particles arriving from one side of the line element (that is, particles 
with X < Xo), then a runs from —n/2 to n/2. 

Py^-"^^ N{\)\d\. (16) 

[We note that the simple arguments found in many physics text books for kinetic theory viscosity in a linear shear flow, 
do not correctly take account of the fact that the particle velocity distribution is isotropic in the frame of the fluid at the 
point of emission (i.e. last collision), rather than at the point where the momentum flux is measured (see for example, Jeans 
1940). In a linear shear flow, this error only changes the coefficient in front of the viscosity coefficient and retains the correct 
form of the equation for the resulting viscous stress. In the circular case, however, we shall find that it is necessary to treat 
this subtlety correctly, even in order to obtain the correct functional form of the viscous stress]. 



2.2 Circular shear flow 

We now carry out essentially the same analysis, but this time for a circular shear flow. We still work in an inertial frame, 
(that is, a frame in which the particle trajectories are straight lines) but in this case the underlying fluid flow is, in cylindrical 
polar coordinates [R, (j)), of the form 

u^{0,m{R)). (17) 

We consider the stress acting on a (small) line element, length dl, centred at point S which in Cartesian coordinates, 
centred at the origin O {i.e. at R = 0), lies at {Ro, 0). The line element lies in the azimuthal direction, and so has unit normal 
in the radial ii— direction. As before, we work in the inerttal frame that co-moves with the mean fluid flow at S and again 
consider the limit that the shear across a mean free path A is much less than the random velocity c {i.e. Rfl'X <^ c).* We 
note that, since RQ,' ~ Cl, it is also the case that A <C -R. We defer until the end of this Section a discussion of the constraints 
placed on A by our neglect of the curvature of particle orbits between collisions. 

As before, we consider particles that are emitted in the vicinity of a point E, located at distance A from S where ES 
makes an angle a with the inward directed radius vector at S. The only difference from the foregoing analysis is that since the 
streamlines are now circular, the velocity of the mean flow at E with respect to that at S has non-zero components in both the 
X— and y— directions (see Figure 2). Specifically, the point E is located at radius -Re = -Ro — A cos a (to 0{X/R)) where the 
velocity of the mean flow (with respect to the origin) is in the azimuthal direction and of magnitude -Rofio — A cos a{Qo + Ro^') 
(Slo and Q' being respectively the angular velocity and its gradient evaluated at S). The radial vector at E makes an angle 

(j) ~ -^sin a (18) 
R 

with the x— axis. Consequently the x— and y— components of the mean flow velocity at E with respect to a frame co-moving 
with the mean flow at S are (to 0(A/-R)) 

^ For a Keplerian accretion disc the disc thickness is -ff ~ c/f2, and thus this approximation impUes that \ <^ H. Thus our analysis 
does not cover astrophysical situations in which the mean free path is comparable to, or larger than the disc thickness (Paczynski, 1978; 
Brahic, 1977) 
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Figure 2. Schematic of circular shear flow. Symbols as in Figure 1. Note that the dashed arrow (the velocity of the mean flow at E 
relative to that at S) now has a component in the x direction due to the geometry of the mean flow 



~ Afio sin Q, (19) 

and 

11^°' ~ -Acosa(no + -Rof^')- (20) 
As before, we can establish the relationship between the angles a and 9 using 
WiT' + c sin 

tana^^^- -, (21) 

Vx + c COS a 

so that, in this case, 

tanS ~ tana + {Qq + Ro^l')\/c + XQq tan^ a/c, (22) 
cos 6 ~ cos a — sin a cos^ a(Qo + Roil')X/c ~ Xflo sin"^ q/c, (23) 
and 

sinO ~ sina + cos^ a(f2o + RqO!)\/c + \Q,q slv? Qcosa/c. (24) 

We note the similarity of these three equations to equations |S| - |H1 since in the circular flow U' corresponds simply to 
(fio + Ro^')- In each equation there is, however, an additional term (proportional to fio) which takes care of the fact that 
the mean flow is not plane-parallel in this case, and that hence each emission point has a non-zero velocity with respect to S 
along the radial vector at S (see equation 119^ . This modiflcation in the relationship between a and 9 is crucial in explaining 
the different form of the viscous stress in the circular case. 

We also need to modify the relationship between L and A compared with what we had previously, which we can do most 
simply by equating the time of flight ES (L/c) with the time (A cos a/(ccos 61 -I- w^"')) required to traverse the distance Acosa 
that separates E and S in the a;— direction. After some algebra we obtain the expression 

L A(l — sin a cos a_Rof^''^/c)~^, (25) 

and hence flnd in this case that 

d/ cos q(1 — sin a cos a_Rof2'A/c) , , 

2^^A • 
Finally, the velocity of each particle emitted from E is 

Vy = csinO + Vy"^ , (27) 

and hence, to the order to which we are working, 

~ csin a — sin^ a cos a 7?o^^'A (28) 
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Figure 3. Velocity components for two particles, 1 and 2, in a linear shear flow originating at emission points located at ±a. Key as in 
Figure 1. 

Comparison of equations II II and 1131 with equations 1261 and 1281 shows that the expressions are identical except in as much 
as U' in the hnear shear flow is replaced by Ro^l' in the circular case. Thus (following the same procedure for finding py as 
previously - that is, by integrating contributions to the y-momentum flux from all points with R < Ro) we obtain (by analogy 
with equation 



We note that this analysis recovers the 'correct' answer (as given by the usual Navier-Stokes argument) that whereas the 
viscous stress in a linear shear flow depends on the velocity gradient, in the case of a circular flow the viscous stress instead 
depends on the gradient of angular velocity. 

Finally, we turn to the issue of our neglect of the curvature of particle orbits between collisions. This neglect is always 
justified, for example, in the case of laboratory Couette fiow, where the fiuid is not subject to external long range forces. In 
the case that the acceleration of the mean flow is provided by a long range central force (i.e. by a force that is experienced 
by the particles between collisions, as in the case of a Keplerian accretion disc) then the change in (x-) velocity of a particle 
between collisions is 5v ~ RQ,^X/c. If (Su ^ c we can repeat the above analysis by adding Sv to v^"^ (equation 19) and hence 
modifying the subsequent equations relating 6 to a. We flnd that (to 0{5v/c)) the relation between L and A (equation (25)) 
and hence the expression for / (equation (26)) is unchanged by this addition but that Vy (equation (28)) now contains an 
additional term ((5« sin"^ a/ccos a). This term however makes a zero contribution to the y- momentum flux, when integrated 
over a, since the contributions from ±a cancel. We thus find that our analysis is independent of the nature of the central 
force providing the acceleration of the mean fiow , provided that Sv <^ c. . This requirement translates into a condition on A 
for a thin Keplerian disc of the form A <g /R, which is more stringent than the condition (A <^ H) derived above in order 
to justify our approach of expanding quantities to first order in Af^o / c. 

3 DISCUSSION 

We have shown that it is possible to obtain the correct expression for the momentum transfer in a both a linear and a circular 
shear fiow using simple kinetic theory. This correct expression implies that an isotropic viscosity transfers momentum down 
an angular velocity gradient. We have noted that it is essential to take proper account of both the geometry and the fluid's 
isotropy. 

Our analysis allows us to understand the flaw in the simple heuristic argument applied to a circular shear flow (described 
in Section 1), which leads to the conclusion that since angular momentum is conserved along particle trajectories, it is the 
gradient in angular momentum which particle mixing tends to smooth out. Particles arrive at the reference patch, at point S, 




(29) 



® Note that it is not necessary that 5v be much less than the shear velocity across a particle mean free path , and that indeed this 
latter condition cannot be met in a thin Keplerian accretion disc. We however find that the symmetry of the terms in the momentum 
flux involving 5v ensures that their contribution can be neglected, compared with those involving Af2o, even when 5v 3> Af2o 
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Figure 4. Velocity components for two particles, 1 and 2, in a circular shear flow originating at emission points located at zba. Key as 
in Figure l.In this example, RQ' + Q = (i.e. the speed of the flow is everywhere constant) but nevertheless the total velocity vector 
(open arrowhead) is larger for particle 1, implying that there is a net flux of y momentum at S. 

with momentum along the streamline (j/— momentum) that derives from two sources - (i) the j/— component of the random 
emission velocity (assumed isotropic in the local rest frame of the emission point) and (ii) the j/— component of the streaming 
velocity of the emission point relative to the reference patch. 

The usual argument given to explain the momentum transfer in a linear shear flow is simply that (ii) has the same sign 
for all the emission points on one side of the reference patch and so that momentum is added to, or subtracted from, the 
reference patch just depending on the sign of the velocity gradient in the mean flow. In fact, there is another contribution, 
equal in magnitude and sign to that described above, which may be understood by considering the particles that are incident 
at the reference patch from ±a. From Figure 3 it may be seen that the a;— velocity (and hence flux) of particles on the +a 
side exceeds that on the —a side {i.e. there is a greater flux of particles at the reference patch whose emission is prograde 
than retrograde). Since the y— momentum from (i) far exceeds that from (ii), this slight asymmetry in the fluxes of prograde 
and retrograde particles produces a net momentum flux that is equal to the more obvious source of momentum flux described 
above. 

We may now apply the same considerations to the circular shear flow, where once again there are the two contributions 
to the 2/— momentum - (i) and (ii) described above. The important difference now is that the mean streaming velocity is no 
longer in the j/— direction. The geometry of the circular arc ensures that this contributes a positive (negative) x— velocity for 
Q > (a < 0), respectively. As may be seen from Figure 4, particles that are incident from +a (prograde particles) have 
a larger amplitude of both Vx and Vy than those from —a. This ensures that the mean j/— velocity of particles arriving at 
the reference patch from a particular emission streamline is not equal to the mean j/— velocity of particles on that emission 
streamline. This disparity is because of the relative boost in the arrival flux of particles whose random velocity is prograde in 
the frame of the emission streamline, compared with those that are retrograde. 

This now allows us to understand the behaviour of a Keplerian accretion disc at a qualitative level. The mean angular 
momentum of particles orbiting at radii less than that of the reference patch is less than that of the reference patch. If the 
distribution of angular momenta of the particles arriving at the reference patch merely reflected the distribution of angular 
momenta of particles on their parent emission streamlines, this would imply that the arrival of particles at the reference patch 
from smaller radii should exert a spindown torque. Instead, the relative boost in the arrival rate of particles that are emitted 
in the prograde direction, ensures that the average angular momentum of the particles arriving at the reference patch exceeds 
the average at the parent streamline. In the Keplerian case, this relative boost in the arrival rate of the prograde particles is 
enough to reverse the sign of angular momentum transfer. Thus, particles arriving at the reference patch from smaller radii 
exert a spin up torque, as required. 

ACKNOWLEDGMENTS 

CJC gratefully acknowledges support from the Leverhulme Trust. JEP gratefully acknowledges continuing support from the 
STScI Visitors' Program. We acknowledge useful discussions on this problem with Peter Scheuer and David Syer. We are 



8 C.J. Clarke, J.E. Pringle 



indebted to the referee, Takuya Matsuda, for interesting comments and for pointing out several typographical errors in the 
manuscript. 

REFERENCES 

Batchelor, G.K, 1967, An Introduction to Fluid Dynamics, Cambridge University Press. 

Brahic, A., 1977, A&A, 54, 895 

Bretherton, F.P., Turner, J.S., 1968, J. Fluid. Mech, 32, 449 
Hayashi, E., Matsuda, T., 2001, Prog. Theor. Phys., 105, 531 

Jeans, J., 1940. An Introduction to the Kinetic Theory of Gases, Cambridge University Press, p. 157 
Kumar, P., Narayan, R., Loeb, A., 1995, ApJ.,453, 480 
Lynden-Bell, D., Pringle, J.E., 1974, MNRAS, 168, 603 

Madej, J. Paczynski, B., 1977 in 'The Interactions of Variable Stars with their Environment', Proc. lAU Coll. 42, eds. R. 
Kippenhahn, J. Rahe & W. Strohmeier, Bamberg: Remeis-Sternwarte, p. 313 
Paczynski, B., 1978, Acta Astr., 28, 253 
Pringle, J.E., 1981, ARAA, 19, 137 

Torkelsson, U., Ogilvie, G.I., Brandenburg, A., Pringle, J.E., Nordland, A., Stein, R. F., 2000, MNRAS, 318, 47 



